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ABSTRACT 

For lensing galaxies we introduce a formal velocity dispersion a^^^^, based on enclosed 
mass and the virial theorem. This is calculated from an ensemble of pixelated lens 
models, and found to be fairly model independent. A sample of 18 well-known early- 
type lensing galaxies and two clusters is found to be consistent with (Ji^^s = fobs- 
Both the early-type lensing galaxies and the clusters can thus be determined as being 
virialized. In a second step we calculate the I-band luminosity and the total mass 
content for the sample of lensing galaxies, which enables us to analyze the mass-to- 
light relation L oc M". We determine a = 0.70 ± 0.08, excluding constant M/ L and 
consistent with previous studies of the Fundamental Plane. Additionally we verify that 
this relation does not extrapolate to clusters, which have a much higher M / L. The 
sample used for this analysis comprises 9 lensing galaxies from the Sloan Lens ACS 
Survey (SLACS) and another 9 from the CfA- Arizona Space Telescope LEns Survey 
(CASTLES) as weU as the lensing clusters ACQ 1689 and ACQ 2667. 

Key words: gravitational lensing - galaxies: fundamental parameters - galaxies: 
elliptical and lenticular, cD - galaxies: clusters: general 



1 INTRODUCTION 

The Fundamental Plane (FP) for early type galaxies is a 
well-known scaling relation between the effective radius Rbh, 
the kinemat ic velocity dispersion ggba and the surface bright - 
ness I<R^ff (|Diorgovski fc Davislll987l : iDressler et al.lll987l ). 
Not so well understood is the mismatch between theoretical 
predictions for the FP on one hand and observations on the 
other. Combining the virial theorem 



formulas. 

With gravitational lensing as an independent measure 
of mass one can suitably analyse the structure of the FP 
as already done in different approaches. iBolton et al.l (l2007fl 
linked l ensing mass and v irial mass, whereas Rusin et al.l 
(|2003al b lTVeu et al.l (l2006ll,|jiang fc Kochanekl l|2007l ) fhere- 
after [JK07]) and lPerreras et al.l ( 20081 ) analvzed the mass- 
to-light dependence 



M oc -RotfO" 



(1) 



M" oc L, 



(4) 



. and the universality of light profiles L oc -Roff^ while assum- 
ing a constant mass-to-light ratio yields 



r> 2 r-l 



(2) 



which we name the Vanilla Plane. In contrast observations 
show a relation with slightly different power indices a and 
6, as in 



which is a representation of the FP. 

By repeating the step from Equation (lU to (0 for the 
more general definition of the mass-to- light relation in Equa- 
tion Q, a and b in Equation ((3)l can be expressed in terms 
of the power index q, which is what we need to compare the 
results with previous FP type studies (listed in Table [B2)) . 
Equating the now a-dependent exponents of a and / yields 



-Rcff oc ct"/'' 



(3) 



a{b) = -2(1 + 26), 



(5) 



with a 1.2 and b ^ -0.8 (e.g.. |j0rgensen et aLlfiggel ). The 
power indices are thus not in agreement with the Vanilla 
Plane indices (a, 6) — (2,-1) for constant M/L, which is 
suggestive of an underlying regularity beyond the above 
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which only applies for a not unique mapping from (a, 6) to 
Q assuming Equation Q. 

In this study we combine both the virial approach and 
considerations including luminosities by means of lensing 
masses Mio„a from 18 early type lensing galaxies and 2 clus- 
ters discussed in detail in Section [21 
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An important role is played by a formal velocity disper- 
sion, which we define as 

For an isothermal sphere this is exactly equal to a line-of- 
sight velocity dispersion. A short introduction to the deter- 
mination of (7iens(^) with the mass reconstruction method 
PixeLens and a detailed description of the lensing sample 
(see Table IBip will be presented in Section [2] We define 
Miona and Mvi, by comparison with Equation Q as: 

M,„. = |gi?„«aL. (7) 

and Mvi, = |g iJettCTob, . 

This in hand we consider the following questions as a 
rephrased puzzle of FP: 

(1) Is the lensing inferred velocity dispersion (Tions from 
nonparametric mass reconstruction equal to the kinematic 
velocity dispersion (Tobs? 

(2) Is this applicable to cluster scale lensing objects? 

(3) Are the computed Mio„s and Afvir consistent with the 
FP? 

(4) Does the FP relation extend to clusters? 

Bearing this in mind we want to give a short overview 

of previous findings. 

In the above mentioned approach bv lBohon et all (j2007l ) to 
the FP problem, they find that ffiena ~ gobs without ta k- 
ing advantage of luminosities (see also lBolton et all lj2008h l. 
This result is comparable to the findings prese n ted in Sec- 
tion |3] of this paper. Furthermore iBolton et al.l |20o3) take 
no baryonic information into account, but a different plane 
is introduced, which emerges from a dimensional change in 
the FP space from surface brightness / to surface density E, 
giving 

-Roft OC cr"™ E*""' . 

But in fact this scaling relation, named a more fundamen- 
tal or mass plane (MP), can be transformed into the shape 
^^^-6771^^-^(1+26™) ^ ^iiicii is consistent with our the- 
oretical assumptions of Equation |6] and thus represents the 
virial theorem. We like to point out that the change to E in- 
troduces basically a redshift dependence, which comes along 
with a grave selection effect that reduces the significance of 
the scaling relation. Moreover this relation is compared to 
the existing FP by introducing a new parameterization of 
the lensing mass 

Ml =c(G~VL(i?„ff/2))', 

where c denotes a str ucture constant an d 5 a newly 
introduced power index. iBolton et all (|2007l ) find that by 
doing so "the tilt relative to the virial relation is essentially 
eliminated" . But basically c and S are again consistent with 
Equation |6] Upon choosing i5 = 1, c becomes lo g37r/2 and 
the new parameterization in l|Bolton et al.|[2007l ) turns into 
a test of the virial theorem. 

Thus the decreased scatter for a MP is rather a natural 
consequence of the added fitting parameter and selection ef- 



fects than a more fundamental scaling relation. Implications 
on structure variations are hardly possible. An appropriate 
treatment on the search for reasons for a tilt in the FP 
originating in certain structural peculiarities includes more 
elaborate approaches that allow for a distinction between 
for instance anisotropy and mass-dynamical structure. 

iRusin et al.l (|2003al ) introduced a self-similar mass 
model for early- type galaxies, consisting of two compo- 
nents: a concentrated component, which traces the light 
distribution and a more extended power-law component, 
which represents the dark matter. They found a strong 
dominance and therefore used the velocity dispersion aiso 
for an isothermal model as a surrogate in the FP yielding 

a mass-to-light relation of Afj^^^, oc L, which was the 

first such result from strong lensing. The error of the slope 
already excluded a constant M/L. While substituting aiso 
they are effectively assuming the virial theorem. 

[JK07] constrain the average stellar mass fraction of 
a halo in favor of adiabatically compressed halo models 
by taking a sample of early type galaxies which consists 
partially of lensing galaxies used in this sample. By means 
of a two-component model stellar and virial mass are fitted 
separately and an isothermal density profile is assumed. 
The paper takes advantage of already K-corrected B-band 
magnitudes and lensing masses and is, because of its com- 
mon subset of gravitational lenses, directly comparable with 
our data. Although it is not explicitly calculated in their 
paper by taking their data we found Mtf^"-^^ oc L, which 
is in p erfect agreement with the result from iRusin et al.l 
l|2003al ). 

Switching from the lensing point of view considering 
Miens to the observational one considering Mvi, enables 
one to compare the FP from previous studies which were 
inferred from len sing with the FP based on stellar dynamics. 
iTreu et al. I (|2006i ) analyze the FP by means of virial mass 
and find that the velocity dispersions for their SDSS lens 
sample are well approximated by (Jiso, which holds also for 
our mixed CASTLES/SLACS sample. 

Ther e are studies, e .g. from I Graham fc Collesd (|l997l ) 
and Trui illo et all l|2004l '). which raised hope that a solution 
for the FP tilt is at least partially given by broken structural 
homology leading to strong correlations between Sersic 
index n and photometri c-independent galaxy properties. 
iGraham fc Collesd l|l997l ) fit R^^" profiles and make use 
of the spatial velocity dispersion at spatial effective radius 
to show the i n fluenc e of structural non-homology, whereas 
iTruiillo et al.l (|2004l ) quantify the contribution to the tilt 
caused by variations of n for a wide range of B-band selected 
early- type galaxies. The results always show that taking 
account of non-homology shifts the FP parameters closer 
but never fully matches the virial expectations. However by 
comparing the M/L — a relation of 25 E/SO galaxies from 
the SAURON sample w ith predictions and virial estimates 
ICappellari et al.] |200^ find that the FP tilt is exclusively 
due to a real M/L-variation, while structural and orbital 
non-homology has a negligible effect, a result also verified 
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in this study. 

Furthermore, progress in estimating Mtot/Mstei was 
recently made by comparing stellar population models with 
the nonparametric mass profiles also used in this paper, 
which allow for scanning the dark matter distribution 
within a galaxy (jPerreras et alj |2005| . |2008| ). They found 
that low-mass galaxies have only little dark matter content 
at all observed radii. On the contrary high-mass galaxies 
have little or no dark matter inside the effective radius but 
at large radii they are clearly dark matter dominated. No 
kinematic and virial assumptions were required. 

In this paper we are using combination of kinematic, 
photometric and lensing inferred data to answer the afore- 
mentioned puzzle. The comparison of (Jobs and aions, which 
is proportional to a comparison between virial mass and lens 
mass, is adequate for answering the above questions 1 and 2 
as will be shown in Section |3] In Section |4] we compute the 
luminosities of the lenses, Mions and Mvir needed for L oc M" 
and check the consistency with other FP studies in particu- 
lar with data from [JK07], who use a common subset of lens 
systems. Subsequently the a-b-parameter plane is generated 
including a wide range of recent FP studies. In Section[5]the 
conclusions are presented. 



2 LENSES AND LENS MODELS 

In this section we introduce the lensing sample and subse- 
quently discuss the lens modelling. 

Our sample consists of 9 lensing galaxies from SLAC^ data, 
9 from CASTLE^ and 2 lensing clusters. We select these 
galaxy lenses using two criteria: 

(1) the lensed images were either point-like sources or 
contain nearly point-like features, and 

(2) the availability of (Toba data. 

Two cluster lenses with such properties are also 
included for comparison and contrast, since previous to 
this paper FP studies were carried out for small scale an d 
large scale objects combined (e.g., ISchaeffer et al.l Il993h . 
As an additio nal motivation it is worth mentioning that 
IZwickvl (|l937l ) originally introduced gravitational lensing as 
a method to estimate masses of galaxy clusters. 

The 9 CASTLES lenses turn out to be a relatively 
inhomogeneous sample, a consequence of the fact that they 
spread over a large range in redshift and effective radii as 
well as lens radii. 

The doubly imaged systems among the CAS- 
TLES lenses are CFRS03.1077, HST15433 and 
MG2016+112. The effective radius of CFRS03.1077 
is not known; hence it is used for the analysis in Section 
[3] but not Section ID In IIST15433 there is a neighboring 
galaxy, but this is thought to only modestly perturb the 



www.slacs.org - The full set includes about 70 lenses, but image 
data was only made available for a small subset. 
^ cfa-www.harvard.edu/glensdata/ 



estimated mass, according to [JK07]. 

Q0957+561 is a special case, as there is a doubly 
imaged galaxy component in addition to the famous 
double quasar. The lensing galaxy is part of a cluster 
that contribut es significantly to the large image separation 
(|Garrett et al.l 11992 ). Consequently a position below the 
general trend in a mass-to-light analysis is expected. This 
lens is an excellent example for the consequences of possibly 
yet unknown image systems. The considerations following 
in Section |3] for (Tic„a-crob3 are carried out for two different 
image configurations: on one hand the 2 double image 
system shown in Figure [1] and on the other hand a single 
double system. 

The quads from the CASTLES catalogue arc B0047- 
2808, PG1115+080, HST14176, B1608+656 and 
Q2237+030. B0047-2808 appears to have a double com- 
ponent source but this is probably not important for macro 
models. PG1115+080 has measured time delays and it is 
also part of a group, which contributes with a significant 
external shear. Neither of these were used for the models 
of this paper; if these are included the lens models tend 
to become rounder, but (Jic„a changes only by 5% to 10%, 
which is insignificant for the present study. For the complex 
lens B1608-I-656 the measured time delays are used and 
make a more significant difference to the lens models. The 
lens Q2237-I-030 is actually the bulge of a barred galaxy. 
In the present study the bulge is treated as an early type 
galaxy. HST14176 is a part of a cluster, which is not 
included in the models. Another possible problem is a large 
uncertainty in the effective radius. 

Several of these lenses have been studied individually 
in great detail. Different paper s sometimes disagree on 
the slope of the profile (e.g., iTreu fc Koopmang I2OO2I : 
iRead et all |2007| , for PG1115-I-080), but agree on the 
enclosed mass. Hence the effect on crions would be small. 

SLACS lenses populate a redshift range from 0.05 to 
0.5. Due to smaller mean effective radii and lens redshifts, 
as a consequence of a limited aperture (3") radius of 
SDSS fibers, the sample of 9 SLACS lenses appears to be 
more clustered in mass-to-light plots than the CASTLES 
sample and has therefore a smaller RMS. The SLACS 
lenses we use are a subset of the ful l SLACS sample fo r 
which point-like features are identified (jPerreras et al.|[2008l l . 

The doubly imaged systems among the SLACS 
lenses are J0037-094, J0912+002, J1330-014 and 
J2303+142. All of them are quite typical with biases 
according to their observation method. J1330-014 is the 
nearest lens with z — 0.08 and shows the smallest (Tobs. 
Among these lenses J0912-I-002 takes a special position. 
It consists of two long arcs which are represented in this 
work as four doubles. Moreover, this lensing galaxy has the 
highest (Tobs among all lenses of our sample. 

Quadruply imaged systems are found for the SLACS 
lenses J1205+491, J1636+470 and J2300+002. Apart 
from J1636-I-470 all SLACS lenses have a larger kinematic 
velocity dispersion than lens velocity dispersion. The mean 
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kinematic velocity dispersion (aoba) of our SLAGS sample 
is 10% higher than {a^^,) of CASTLES lenses. 

Also J0737+321 and J0956+510 are thought to 
be quads, but in each of these only 3 images are used, as 
the astrometry of the faintest image was too uncertain. 
J0737+321 is with z = 0.32 the most distant in our SLAGS 
subsample and belongs to higher z lenses in the whole 
catalogue. 

Finally we consider the two lensing clusters AGO 
1689 and AGO 2667. 

AGO 1689 has a very large number of multiply imaged 
systems found by Broadhurst ct al. (2005). In the present 
work this cluster is modeled by a set of two 5-image sys- 
tems, six 3-image systems and one doub le. The additiona l 
systems are known to affect only details ()Saha et al.ll2007h . 
Note that there are many more imaged sources, but adding 
those to the model does not change Mi„„b, i.e. the mass 
model is tightly constrained by this set of image systems. 
The kinematic line-of-sight velocity dispersion cTobs = 1400 
km s~^ of galaxi es within the cluster was taken from 
iLokas et al.l ()2006l ) for a subset of 130 galaxies in the 
inner region of the cluster with velocities \v\ < 3000 km 
s~^, which contains most likely the biggest mass fraction 
responsible for the lensed images. This average value 
applies for a radius of around 400 kpc, a region where the 
formal velocity dispersion seems to be sufficiently flat and 
in which roughly hal f of the pro j ected radii of the 130 
galaxies considered in iLokas et al.1 l|2006h are to be found. 
Furthermore the value is not too far away from the Einstein 
radius or outermost image position of around 240 kpc. In 
order to estimate the I-band magnitude of the cluster, the 
130 brig htest out of 840 galax ies are taken from a cluster 
survey of lMolinari et all (|l996t ) for which the Gunn g,r and 
i magnitudes were provided. Together with K-correction, 
evolution correction and galactic extinction we obtain 
Li = 2.82 X W^^Lq. 

AGO 2667: For this cluster three 3-image systems and 
one double were known to derive the formal mass-only 
related velocity dispersion curve. The kinematic velocity 
dispersio n cTobE, = 960lj!^n km s~^ of this lens was deter- 
mined bv lGovone et al.1 (j2006h from a sample of 21 galaxies 
in the inner region of the lensing cluster with a radius of 
IWh^Q kpc, which is in the same order of magnitude as 
Riaas = 98h^2 kpc. However, since photometric data for 
estimating the total flux of galaxies within the cluster was 
not available, AGO 2667 is not included in the mass-to-light 
plots of Section |4] and consequently there was no need to 
determine -Rett for the mass estimate. 

The a bove lensing data are mo d eled using th e PixeLens 
progranjfl (jSaha fc WiUiamsl l2004l : iGolesI l2008l ). PixeLens 
reconstructs the projected mass in a pixelated manner by 
solving a set of linear constraint lensing equations on the 
mass distribution by means of the given image positions, 
the redshifts of lens and source, the Hubble time (herein 
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E 

= 100 - 



10 20 30 40 

radius (kpc) 

Figure 1. Left panel: projected mass distribution of the CAS- 
TLES lens Q0957-I-561. The box size is 4 X 4 arcsec^. The dots 
mark two doubly imaged systems. Right panel: formal velocity 
dispersion crions. The vertical dashes mark the radial position of 
the lensed images. The same curve can be seen as second from 
top in Fig. [2] 



/i = 0.72 is always assumed) and optionally the time delays 
between the lensed images. There are certain requirements 
for the mass distribution to be fulfilled. It has to be 
non-negative, centrally concentrated, with a local density 
gradient pointing less than 45° away from the center, 
inversion symmetric (optional), it must not have a pixel 
which is twice the sum of its neighbors except possibly the 
central pixel and the circularly averaged mass profile needs 
to be steeper than i?""'^, where R is the projected radius. 
Hence an underdetermined set of equalities and inequalities 
is obtained. Subsequently a Monte-Garlo approach is used 
to sample over the mass map in order to determine an 
ensemble of lens models, for which the ensemble-average 
appears to be the best single model representation. Of 
course all the uncertainties on any parameter can be derived 
from this model ensemble. 

Since PixeLens has been extensively tested in other pa- 
pers, we won't get into details, but two points are worth 
noting: 

( 1) For tests of th e recovery of simulated galaxy lenses, 
see lRead etahl (|2007h . 

(2) For the recovery of gross features of even extended 
lens structures from the informa tion encoded in the im - 
age positions o f lense d objects, see I S aha fc Wfihamsl l|200ll ): 
iFerreras et al] l|2008l l. 

For each lens an ensemble of 100 mass maps with 
21 X 21 pixels each has been computed, from which the 
mass-profile and therewith the formal velocity dispersion 
(Jicna is derived with an 90% uncertainty, as one can see 
e.g. for the lens Q0957-I-561 in Figure [T] In Appendix iBl we 
compare the average of an ensemble containing 100 models 
with larger ensembles containing up to 10000 models to find 
that already small ensembles are sufficient to determine 
a fairly exact velocity dispersion. Since o-ions(-R) is not 
sensitive to ensemble enlargement the number of models is 
fixed to 100 throughout this analysis. 

Two points shall be emphasized here. Firstly, the error 
bars in the right panel of Figure [T] represent the model 
dependence for an ensemble of 100 models and, as one can 
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see, it is not large. Secondly, it is sometimes stated that 
the enclosed mass M(< R) is known for R = i?Ein, the 
Einstein radius and unknown for any other R, but this is 
oversimplified. In fact M{< R) has some model dependence 
at all R, but is minimal at -REin- c"ions(J? 7^ -RBin) has 
a larger uncertainty than cricna(i?Ein), but is still fairly 
well constrained, as one can see in Figure [1] The velocity 
dispersion at the radial position of the outermost image 
(Tiens(J2iens) as a quautity, which is as well constrained as 
i?Eiii, is basic to this paper. 

The PixeLens input files, mass maps as well as the 
formal velocity dispersion curves can be found in Appendix 
[B]of the online version of this paper. Note that for lensing 
clusters the velocity dispersion of the galaxies on their orbit 
around the center of the cluster is considered instead of the 
stellar velocity dispersion as in the case of lensing galaxies. 
The references for all the lenses and further details can be 
seen m Table [BT] 



- 500 




3 cr,e„s VERSUS (Jobs 

The formal velocity dispersion curves (7ions(^) illustrated in 
Figure [1] for the lens Q0957+561 are now computed for all 
the lenses. Figure [5] shows all these curves except for the 
cluster, which are excluded for the sake of readability and 
their comparatively high (7ic„b values. 

In the following we concentrate on the formal velocity 
dispersion at a radius of the outermost image position 
o"iona(i?iona) aud at thc effectivB radius aians{Rait)- Concern- 
ing the latter we cannot take for granted that the velocity 
dispersion curve at effective radii is still sufficiently flat. 
Because of this, when considering cric„a(i?ott) we exclude 
the lenses Q2237, HST15433, J0737, J0912, CFRS03 and 
J0956, for which this condition is not fulfilled. 

In terms of absolute values the curves for CASTLES 
lenses extend in average to larger radii whereas the curves 
of SLACS lenses are smaller due to a limited aperture of the 
SDSS fibers. Note that 7 SLACS lenses and 3 CASTLES 
lenses show a clear cuspy shape of the formal velocity 
dispersion curve towards inner radii as it is the case for 
the majority of early-type galaxies also in oth er velocity 
dispersion field studies (e.g.. ICoccato et al.ll2009l ). However, 
in some cases anomalous galaxies exhibit a rising velocity 
dispersion profile, whi ch might be re l ated t o the presence 
of a disk according to ICoccato et all (|2009l ). Additionally, 
the pixelated approach causes a variety of diff'erently 
shaped velocity dispersion profiles differing especially in 
central regions. This leads consequently to large error bars 
and a decreased sensitivity in the center, rendering an 
interpretation of the profiles at smaller radii rather difficult. 

It should be emphasized that the comparison of either 

Ciena (ftena) Or (Tlena (i?eft) wlth ffoba mCaSUred wlthlU aU 

aperture, is a proper procedure, since the Ri^hb is in average 
less than a factor of 2 different from the aperture radius 
and for most lenses aions remains unchanged. For J0737, 
J1205, J1330 and J2300 the formal velocity dispersion 



Figure 2. Formal velocity dispersion curves for all the galaxies. 
The grey solid lines denote SLACS lenses, the black dashed lines 
CASTLES lenses. For the sake of readability the error bars, in- 
dicating the range of ensemble models, are shown only for one 
velocity dispersion curve. The radial scale is normalized to the 
radius of the respective outermost image Ricna which is indicated 
by the horizontal solid line. Open circles denote the effective ra- 
dius provided it is located in a fairly flat region of the velocity 
dispersion curve. Q0957 is shown twice in this plot: Q0957-1D 
models one doubly imaged source while Q0957-2D models two 
doubly imaged sources. 

curve ends before reaching the radius of 3 arcsec, that is, 
the mass does not contribute to the lensing effect, but 
nevertheless tJoba can be taken as an indicator for the 
real velocity dispersion. With other words the velocity 
dispersion measurements at aperture radius are probably 
not representetive since the main mass of the lens is smaller. 

Comparing the curves labeled Q0957-1D and Q0957-2D 
shows the probable effect of adding formerly undiscovered 
image systems. As for Q0957 (Tic„a(i?) varies consider- 
ably when a formerly unseen doubly imaged system is 
added. This also affects the relation between (Tion5(-Rion5) or 
Ciena (i?cff) and the kinematic velocity dispersion cToba. 

Both (TicnB(i?icnB) and criona(Jiotf) plotted against iToba 
can be seen in Figure |3l The comparison between the 
observed kinematic velocity dispersions and the mass-only 
related velocity dispersions reveals how virialized the lenses 
are, because aiona ~ Coba is another representation of the 
virial theorem in Equation [S] We constrain the fit by fixing 
it to the (0, 0)-point, because a bias would have no physical 
relevance. For uicna-values at effective radius instead of the 
radius of the outermost image position i?icna the scatter 
around the best fit decreases considerably. Although all 
o"iona(J?otf) are within the error bars of o"iena(^iiena), changing 
the radii for the determination of the criena-croba-relation 
might consequently be the right thing to do, since therewith 
the relation is build on a common basis. 
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Figure 3. trions-f obs-plot for all the galaxy lenses. The filled cir- 
cles refer to formel velocity dispersions cricna measured at outer- 
most image i?icna. The open circles show o'io„s{i?ott). Note that as 
in Figure|2]Q0957 is shown twice for different image systems. The 
dashed (dotted) line represents the fit for the solid (open) circles 
including Q0957-1D (Q0957-2D). 

Furthermore we included the two-double (2D) and the 
one-double (ID) system of lens Q0957 in Figure|3]to demon- 
strate the grave difference in crions of the former outlier, re- 
ducing the RMS in the o-ic„B-o-obB-plot from 55 km s^^ for 
the cric„a(J?ic„a)|i_D-fit to 43 km s""*^ for the (Jic„a(7?ctf)i2D-fit. 
We conclude that generally a more complete lens system is 
to be favored and henceforth we only consider Q0957-2D. 
The linear best fits fixed to the origin for (Tions (i?cft) and 

O"lcn3(i?lona) yicld 

(7ob3 = (1.03 ± 0.05) X a,„„3(-Roff), (8) 
fTobs = (1.04 ± 0.04) X a,„3(-Ric„3). (9) 

As an aside the y error bars p lotted in Figure [3] are th e 
observational error s taken f rom iKoopmans fc Treul l|2003l); 
iTreu fc Koopmand (|2004|): iTonry fc Franxl (Il999|) : Tonrvl 
(Il998h: l Ohvama et al.| j2002l '): ^Koopm aiis et all l|2003h : 
iKoopman s & Trcu (20o 3) and [P ohz ct ajj (Il992l ) for the 
CASTLES lenses and iBolton et al.l (j200a ) for SLACS 
lenses. The x error bars represent the statistical errors 
of the formal velocity dispersion for an ensemble of 100 
models of possible mass distributions. Thus the rather 
small error bars can be understood as a relatively model 
independent lensing mass and formal velocity dispersion. 
The errors are taken from a radius closest to J?ic„3 since the 
pixelated approach only allows for discrete steps in radius. 
One could argue about the significance of these errors, 
because changes in the image positions or lost information 
like additional image systems or mass contamination of the 
light path can lead to fairly different results. 

However, the fits for aabeio'icTis) (Equations (|8I9|) ') make 
clear that a one-to-one correlation between M-i^^s and A/vir 
of the lensing galaxy is probable. It is important to know 



Figure 4. Plot of Ricas/Rdf against CTions/fobH. The lines indicate 
extreme scenarios of formal fits for our sample of early-type galax- 
ies which show large errors. For the (Tic„3(iJctf)-fit of the whole 
sample the positive trend is insignificant. The trend inverts when 
excluding MG2016. In other words there is neither correlation nor 
anti-correlation, meaning that in average the density profile for 
all lenses is consistent with an isothermal ellipsoid. 



whether our sample is dominated by a certain kind of 
model far from p{r) ~ corresponding to a constant 
(Jicns. For that we can study the correlation between the 
ratios aio„s/o"obs and i?ic„s/-Rctf. In consideration of the virial 
theorem one can state: 

If there is an (anti-)correlation between CTicna/coba and 
Riona/Raii thc dcusity profile p(r) of the lens should be 
(flatter) steeper than . 

Figure |4] shows this relation for both (Jie„B(J?ie„B) and 
f"ions(J?oft). As for the first, the best fit shows a positive 
trend with large error bars. For (Jic„s(J?oft) the positive trend 
is insignificant and the opposite result is not excluded by 
the error bars. By neglecting the outlier MG2016 with a 
possibly underestimated R^i, as will be discussed in Section 
|4l one finds the inverse trend to be likewise significant. 
Nevertheless it needs to be emphasized that by excluding 
only one of the labeled outliers in Figure |4] the slope is 
strongly affected and can change its algebraic sign. Thus 
we cannot retrieve a strongly significant statement. In 
such exclusion scenarios we obtain slopes consistent with 
constant cr-ratio. Our sample of early type lensing galaxies 
for (Ji;,„s(i?i;,„s) ((Tic„a(i?etf)) is clustering around a mean of 
1.3 ± 0.3 (1.6 ± 0.3) in 7?ie„,/-Rctt and around 0.91 ± 0.04 
(0.96 ± 0.06) in o-ie„,/cr„b, excluding MG2016 because of 
its extraordinarily high 7?io„B/7?Bft-ratio. Since we cannot 
find any type of correlation throughout our sample we can 
summarize that ctio„s is model independent. 

Extending the (Tobs-o"icns-plot in Figure [3] to 100 kpc 
scale, as to be seen in Figure [S] we can find that lensing 
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Figure 5. Like Figure |3] but with the two clusters AGO 1689 
and AGO 2667 included. The straight hne fits do not include the 
clusters. The relation between (Toba and cicnu extends to cluster 
scales. 



clusters fit quite well to tlie previously found correlations 
for aobs{o'icas) (Equations (|8l9p ). 

Going from an assumed isothermal profile to a 
Hernquist profile makes the velocity dispersion of Equation 
((SJ change to sho wn in Equation (lAlfl of Appendix 
{K\ along the lines of iHernquistI (|l990l ). This step yields 
a change of less than 19% of (Tic„b for most lenses, apart 
from few exceptions like P1115, which turned out to be 
an outlier already in Figure 3] Furthermore cluster AGO 
2667 shows a velocity dispersion increased by 33%. In 
general central regions o f galaxy c l uster s are best fitted 
by a Hernquist profile l|Hernauistl Il990l ) for the stellar 
component of the in ner cD galaxy and a NFW model 
avarro et al.l []"9 96'l for the dark m atter component, 
shown by e.g. iPadmanabhan et al.l (|2004 ). That is why we 
can expect significant changes going from (Tic„3 to on 
larger scales. However, fitting the (Ticna-cToba relation for a 
Hernquist profile as done before with an isothermal model 
for (Tiona reveals a slightly steepened slope compared to 
Equation © of (1.13 ±0.04). The clusters still agree to this 
relation within the error bars. 

The dynamical state of galaxy clusters is hard to deter- 
mine. There are many contradictory investigations on this 
topic. Optical and X-ray data on one hand indicate ongo- 
ing formation processes on substructure level (e.g., ISteinI 
Il997l : ISolanes et ahll 19991 ) , which should be considered in es- 
timates of Mvir. On the other hand statistical comparisons 
of different mass estimates from optical and X-ray obser- 
vations and weak lensing show perfect agreem ent on scales 
much greater than the core radius 7?core (e.g., IWu fc Fand 
1 19971 '). Still, on scales of core radii there are discrepancies 
between X-ray and mass measurements by means of weak 
lensing. lAUen l|l998l ) suggests to consider substructure and 



line-of-sight alignments of material towards the cluster cores 
since they will increase the lensing masses without affecting 
X-ray data and to take account of the dynamical activity 
which mi g ht cau se the X-ray analyses to overestimate -Rcoro . 
IXu et all (|2000l ) take this apparent dichotomy as an indica- 
tor of the transition from previrialization to virialization. In 
this paper however we can probe the virialization state for 
the two clusters at -Riona, which is in both cases not far away 
from -Rcorc- The core radii of the X-ray selected AGO 2667 
and AGO 1689 ar e abou t (76 ± 8) kpc (|Govone et al.ll2006l ) 
and (80± 15) kpc lAllenI (l998) respectively. Thus with (Ji„„a 
at i?otf ~ 98 kpc for AGO 2667 we already probe the core 
region. For AGO 1689 R^h is roughly 238 kpc, which is 3 
times the given core radius. By adjusting to smaller scales 
Ciena (J^coro) bccomes ~ 1000 km s~^ and marginally fails the 
relations (O and It should be emphasized that unlike 
the sample of lensing galaxies AGO 1689 iicore is not in a 
sufficiently fiat region of (Tic„a and thus not comparable with 
the relations for which this was a requirement. Since strong 
lensing unveils mass regardless of underlying dynamics one 
can summarize that also in view of findings from previous 
studies clusters in a wide range of radii can be regarded as 
virialized. 

Nevertheless, the correlation between the kinematic ve- 
locity dispersion (Toba and (Jic„a is hard to decipher. First the 
scatter around a best fit that is smaller (larger) than the 
scatter around the FP in the {R^tt,a,I) parameter space can 
be understood as a hint on a basically mass dependent (stel- 
lar dynamics dependent) (Toba. Of course it can also be seen 
as a merely statistical scatter that is infiuenced by a possi- 
bly biased lens sample. This allows for drawing the following 
conclusions: 

(1) The small scatter and the slope of the best fit of ~ 1 
makes criona a good surrogate for (Joba, which is independent 
of a particular density profile model. 

(2) The included elliptical galaxies are thus virialized and 

(3) the relation can be extended to larger scale objects 
like clusters, as we can see in Figure [S] 

With this in hand we now want to analyze the mass- 
to-light relationship for the given sample and compare it to 
the governing FP of early type galaxies. 



4 MASS-TO-LIGHT RATIO AND THE 
FUNDAMENTAL PLANE 

As a first step to a mass-to-light relation for this sample of 
early type lensing galaxies we K-correct given I-band mag- 
nitudes (centered on 814 nm) and SDSS-i-band magnitudes 
(centered on 753 nm) in AB units to rest frame I-band 
since they provide the most complete set of magnitudes 
for our sample. These are taken from the GASTLE Survey 
homepag^fl and iBolton et all l|2006l '). In the case of the 
galaxy cluster AGO 1689 we obtain the overall magnitude 
by summing over the fiuxes of the galaxy content using the 
catalogue of iMolinari et al" I l|l996l ). Hence the K-correction 
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Figure 6. Visualization of K Correction: The black solid 
curve shows the flux template of an elliptical galaxy. The grey 
solid line represents the HST WFPC2 I-bandpass taken from 
www-int.stsci.edu/instruments/wfpc2/Wfpc2_thru/. The dashed 
curves are showing the denominator of the integrand in Equa- 
tion (TDll for 3 lenses: MG2016 (z = 1.01, dashed line), B0047 
{z = 0.485, dotted line) and Q2237 {z = 0.04, dash-dotted line). 



is based on SDSS, HST and ESO spectral template^ 

We carry out the K-correction from first-principles 
in preference to a black box program. Following 
lOke fc San dag3 l|l968t) we compute the K-corrected flux ac- 



cording to 



2.51og(l + z 
+2.5 log 



/„°°F(Ao)S.(A) d\ 



/;°F(Ao/(l + ^))5.(A) dX 



(10) 



where Kx denotes the K-correction for the x-Band expressed 
in magnitudes. The band width is smaller in the redshifted 
galaxy, which leads to the first term in (|10p . A source 
spectrum F(A) is redshifted through fixed spectral-response 
bands Sx or bandpasses respectively of the detector. The 
flux at an effective wavelength in the rest frame of a galaxy 
of redshift z, transformed from the effective wavelength Ao 
of the detector by Xo/{l + z), will differ from the flux of a 
galaxy at rest. This leads to the second term in 1101 Figure 
[H] visualizes the denominator of the integrand in Equation 
(|10p . where the I-bandpass is multiplied by the redshifted 
flux te mplate of an elliptical galaxy taken from lKinnev et al] 
l|l996h . As an aside the apparent SDSS magnitudes are on 
an AB basis within 3%, which only leads to minor correc- 
tions and is therefore neglected in the following analysis. 



5 The spectral templates for WFPC2, SDSS ACS 
and ESO telescopes are taken from the STSCI home- 
page www-int . stsci . edu /instruments / wfpc2 /Wfpc2_thru / , 
stsci.edu/hst / acs / analysis / referencc_files / synphot_tables.html 
and filters. Is. eso.org/efs/efsji. htm 



Note that the K-correction is realized with the exact tem- 
plate for A < 570 nm. For higher wavelengths we assumed 
a constant flux for the sake of simplicity. The deviations re- 
sulting from this approximation are even in the worst case of 
a hardly redshifted galaxy in the upper A-range like Q2237 
of only 0.3% for Lj. This leads to negligible corrections for 
all following quantiti es. Furthermore galac tic extinction cor- 
rections according to lSchlegel et al.l (| 19981 ) are applied to the 
fluxes. The luminosities are calculated in units of solar lu- 
minosities according to an AB magnitud^fl /q = 4.57 for 
WFPC2 da ta and = 4.48 for S PSS data calculated along 
the lines of iFukugita et al] l|l995h . Subsequently we correct 
for passive M/L-evolution with a slope of 



= -0.397 



dlog M/Li 
dz 

inferred by stellar popula tion synthesis models taken from 
iBruzual fc CharlotI l|2003h . 

Having the I-band luminosities Li of all lenses in 
units of solar luminosities Lq and the velocity disper- 
sions from Section |3] we can analyze the underlying 
mass-to-light relation. Figure [7] shows the lensing mass 

Aflcna = -RcffcitnaCiilona) aud thC vlrlal maSS Mvir = RalfO'abs 

plotted against I-band luminosity. The plot also provides 
a curve representing a constant M/L or according to 
Equation Q a (a = l)-line respectively. 

A closer look at the V-band luminosities for selected 
galaxies reveals that HST14176, B1608 and MG2016 emerge 
as outliers with mass-to-light ratios < 1. This can be ex- 
plained by nearby groups and clusters (e.g. in the case 
of HST14176) or mass-contamination influencing the path 
of light. Another reason can be uncertainties in the ef- 
fective radii, as already mentioned in Section |3] If we 
take for HST14176 (MG 2016) R^a = 1.06 (0.31) arcsec 
jTreu fc Koopmanj |2004|) instead of the used 0.71 (0.22) 
arcsec (IRusin et al.l l2003bl ) then Miono ~ ^oftcfons would in- 
crease by a factor of ~ 1.5 (1.4), since no grave changes in a 
for a flat formal velocity dispersion curve are expected. This 
leads for HST14176 to a lensing mass of 5.14 x lO" Mq in- 



stead of the former 3.43 x 10^^ Mq which is then also in the 
V-band clearly below the {a — l)-line. If such uncertainties 
are the true cause for comparatively high luminosities, then 
we also need to adjust Mions and Afvir in Figure [7] However, 
changing R^a or excluding the problematic lenses from the 
fit has a negligible impact on the slope a using A'/io„a and 
only small impact using Afvir, changing a from (0.80 ±0.10) 
to (0.84zfc 0.10 ). It should be emph asized that we hold on to 
the dataset of Rusin et alj (|2003lj ). because it provides the 



effective radii computed on a common basis for the whole 
CASTLES subset of our lensing objects. 

Both sets of data points for cToba and for o-ic„a(J?i„„a) are 
fitted for the whole sample and reveal the slopes: 

a = (0.70 ±0.08) for Ar,„„a, 
a = (0.80 ±0.14) for M^i,. 

It shows that a = 1 is in any case clearly excluded. 
Figure [7] shows the best fit for both Afie„a and Afvir- Note 
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Figure 7. Lensing mass and virial mass against I-band luminosity 
for all galaxies. The triangles denote masses calculated using cric„a. 
The squares refer to masses calculated with (Tobo • The best fits for 
oc L are plotted for (Tic„a(i?ic„a) (dashed line) and for iToba 
(dotted line). The solid line refers to a constant A//L-ratio. 



Figure 8. Mass-versus-light plot for data from [JK07]. The black 
circles denote a subset of lenses included in the lens sample of 
this paper. The grey circles are residual lenses. The dotted line 
represents the best fit for the whole dataset taken from [JK07]. 
The dashed line refers to the a-ic„a-fit as seen in Figure [7] 



that the fits in the plot cannot be extrapolated to lower 
masses, which would mean that judging by the intersection 
with the [a = l)-line the luminous mass would overtake 
the total mass content. The plot and therewith also the FP 
of nearby lenses show that more massive galaxies have a 
larger dark matter fraction. 

In Figure [8] the lens sample from [JK07] together 
with a best fit is shown. Their data from stellar-dynamical 
measurements on 22 early-type models contains a common 
subset with the present study. Note that the data in [JK07] 
was given in B-band luminosities, which explains the shift 
of the data points towards lower luminosities in most 
cases. The fit for the whole [JK07] sample yields a slope 
of a — 0.88 ± 0.12. As in Figure [7] a curve indicating a 
constant mass-to-light relation is included. 

From these plots we can summarize, that 

(1) the slope of the best fit for (Toba is consistent with the 
one for the lensing sample of [JK07] within error bars, 

(2) the slope of the fit for uiena is not consistent with the 
fit for data from [JK07], although the error bars do overlap, 

(3) only the fit for the I Jiang fc Kochanekl (|2007l 'l sample 
is consistent with a — 1, 

(4) the slopes of the (Tic„b- and (Toba-fits (for the whole 
dataset and for a reduced or, due to uncertainties in R^h, 
changed dataset) are clearly excluding a = 1 within their 
error bars and thus do not agree with a constant M/L ratio. 

In Figure [U] we extend determined mass-to-light rela- 
tions to larger scales. For cluster size objects like AGO 1689, 
Rett is of course not defined. Nevertheless, one can still 
use the mass quantity Ra^ to compare the mass-to-light 
behavior of early type galaxies and clusters. The kinematic 



line-of-sight velocity dispersion cToba = 1400 km s~^ of 
galaxies within the cluster was taken from iLokas et al.l 
(2006) for a subset of 130 galaxies in the inner region of the 
cluster with velocities \v\ < 3000 km s~^, which contains 
most likely the biggest mass fraction responsible for the 
lensed images. This average value applies for a radius 
of around 400 kpc, a region where the formal velocity 
dispersion seems to be sufficiently flat and in which roughly 
half of the p rojected radii of the 130 galaxies considered in 
iLokas et al.l |2006) are to be found. Furthermore the value 
is not too far away from the outermost image position of 
around 240 kpc. Therewith J?otfcroba(^otf) and /?ctfO"oba(^icna) 
are determined. 

The data points for the cluster deduced from the for- 
mal velocity dispersion ai^as at Ricns and R^it are included 
in Figure (9] As expected neither the relation M" oc L 
with a ~ 0.70 nor with any other slope presented above 
does extend to clusters. As shown bv lSchaeffer et al] (Il993l ) 
galaxy clusters follow indeed a different FP relation. We can 
make up a region in the mass-to-light plot for cluster sized 
objects, which lies far b elow al l previ ous lines and matches 
the findings of Schacff er et al.l (| 19931 ) for a FP consisting 
of 16 clusters. One should keep in mind that for early-type 
galaxies M/L can be a suitable dark matter versus baryon 
estimator because L tracks pretty much all baryons. But 
this is not a good approximation for clusters, whose total 
baryonic mass is generally believed to be made of 80% ho t 
diffuse gas and only 20% galaxies iFukugita et al.l l|l998h . 
In order to correct for this discrepancy one might add the 
missing 80% expressed in terms of luminosity. Hence the 
luminosity of AGO 1689 is shifted to 1.4 x 10^^ Lq. Despite 
of this correction we obtain a value significantly below the 
given fits. Thus clusters can nonetheless be regarded as 
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to b relation and a common a of ~ 0.80, although mea- 
sured separately. This value is verified in this paper by the 
mass-to- light relation found for Afvi,. It can be seen that 
(a, b) for the s l ope o f the aiena-fit i s close to the results of 
iGuzman et all (ll993h:IColless et all l^2OOl^:[J0rgenserl et all 



10^= 10 " 10^^ 10^= 



(Il996t): IScodeggio et all (|l997l ): iLucev et all (Il99ll ) and 
iDressler et al.l (Il987l) in ascending orde r of distance in (a, b) - 
space. Except f or [Hudson et all [l997l): iPahre et all l|l998l) : 
iGibbons etall (|2000^ and iBernardi et al.l (|2003i rtlie errors 
of previous (a, 6)-studies, as far as they were given, overlap 
with the error b a rs in this s tudy. In particular t he results of 
|j0rgensen erall (|l996l ) and lColless et al.1 |20o3) agree with 
the upper limit of a-values from the aiona-fit. However, the 
Q estimate from the da t aset of [JK07], which matches the 
result from iRusin et al.l (|2003bh can be excluded. Since for 
all previous FP type studies kinematic velocity dispersion 
measurements are used, our findings suggest that the real 
underlying (a, b) values are even closer to the lower right 
corner of Figure [TOl 



Figure 9. Like Figure [7] but extended by cluster AGO 1689. As 
expected clusters are not on the FP, because they have a higher 
mass-to-light ratio. Solid symbols denote masses calculated with 
a{Ref{), open symbols denote masses calculated with a{Risns)- 

highly dark matter dominated. 

We can summarize that our results are in good agree- 
ment with most of the recent FP-type studies, as one can 
see in Figure (TO] In the two plots the FP parameter study 
results of the references listed in Table lB2l are presented (left 
panel) together with the results of this paper (right panel). 
Recovering the FP of early-type galaxies by means of the 
photometric-independent (Ji„„a shows that non- homologies 
like structural and orbital anisotropics, which might change 
the photometrically determined central velocity dispersion, 
have sma ll to negligible i mpact on the FP tilt, as also 
shown bv lCappellari et al.l (120061 ). 

The FP parameters of our analysis are determined in 
consideration of the relations a — 2a(2 — a)~^ and b = 
-(2-a)-i: 

a =1.08, 6 = -0.77 for cri™^, 
a = 1.33, b = —0.83 for aoba, 

corresponding to q = 0.70 ± 0.08 and a = 0.80 ± 0.10 
respectively. Upper and lower limit of the (Tie„a-fit are also 
drawn into the plot and exclude plainly the M oc L case of 
the Vanilla Plane. 

Moreover, the FP parameters found in this study are 
conspicuously surrounded by the ones found in other stud- 
ies (see Table IB2[I . For example in recent SDSS results 
for nearly 9000 early-type galaxies in a redshift range of 
0.01 < 2 < 0.3 the parameters ar e determined to a = 
1.49 ± 0.05 and b = -0.75 ± 0.01 (|Bernardi et al.ll2003D . 
and as an aside have no commo n a, since Equation ((S} 
does not hold. On the other hand IDressler et al.l |l98^) in 
one of the first FP parameter studies present parameters, 
which are almost perfectly in agreement with the fixed a 



5 CONCLUSION 

We can summarize the findings of this paper as follows: 

(1) Independent of the details of lens models the lensing 
masses and virial masses basically agree, since (Tions ~ fobs, 
as demonstrated in Section [S] This verifies the virial theo- 
rem. 

(2) The relation between the lensing inferred velocity dis- 
persion (Tiona and the observed kinematic velocity dispersion 
(Toba extends to cluster sized lensing objects within rather 
large uncertainties originating from a poorly defined scale 
radius R^h as shown up for the two galaxy clusters AGO 
1689 and AGO 2667. 

(3) Using the results for (7ic„a (cToba) in Section |4] the lens- 
ing mass (virial mass) is calculated according to Af ~ R^ffa^. 
We find the mass-to-Ught relation M°^J°^°-°^ oc L for the 
whole sample and M^l^''^^'^'^ oc L to be consistent with 
most other FP type studies. We point out that the FP de- 
fined by using cricna(J?icna) is based on lensing velocity dis- 
persions within -Ricna , which is not correlated to the effective 
radius. In order to render the used quantities unequivocal we 
analyze the change in (Ti„„a-(Toba switching from -Ric„a to i?otf 
and find only a marginally different slope, though a reduced 
scatter in the criena-CToba-plot can be seen. A few lenses are 
problematic outliers due to observational uncertainties but 
excluding these does not effectively change the result. With 
R oc (T,^„',Jf J""'^^ the FP of early type galaxies is recovered, 
excluding clearly the Vanilla Plane. Thus also non-homology 
as a reason for the FP tilt can be excluded. 

(4) As shown for AGO 1689, clusters are far from the FP 
since they have a much higher dark matter fraction than 
early type galaxies. 

The FP tilt discovered by IDressler et al.l |l983) and re- 
covered in this study using cric„a as a surrogate is an often 
discussed matter (see Table IB2|I in astrophysics. The rea- 
sons for the deviation from the Vanilla Plane are hard to re- 
solve, because neither the mass-structure, the mass-to-light 
ratio nor the dark matter fraction are directly and indepen- 
dently observable. Until a consensus on the explanation for 
the FP is found it is necessary to focus on quantities which 
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Figure 10. Left panel; the a-b-parameter-space according to Equation JSj. The dashed line represents the mass-to-light power index a 
for related a and b according to Equation l[5)l . The plot shows FP parameter from previous studies referenced in Table IB2I The error 
bars are included as far as provided in the references. Right panel: like left plot but with results from this paper. The squares mark 
upper, lower and mean value of the fit using cicnu for the whole sample of early-type galaxies. The open circle denotes the fit using CToba- 
The black filled circles denote other o-values, like e.g. the Vanilla Plane or the fit for data from [JK07]. For the sake of readability and 
comparison the grey filled circles corresponding to the data shown in the left panel are included. 



are unequivocally related to a certain physical entity. For 
this purpose criens is proposed in this paper, since it fulfils 
the necessary condition of preserving the viriality for both 
elliptical galaxies and clusters. 
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Figure Bl. Absolute deviation of formal velocity dispersion 
curve of a small ensemble from the one of a large ensemble in 
terms of percent plotted against the radius in terms of -Ricna . The 
dashed vertical line denotes i?ic„a . The listed markers indicate the 
deviation of cri^^s at the effective radius. 



Xis) 



for sC s < 1 
for 1 ^ s oo 



according to iHernauistI l|l990l ). 



APPENDIX B: 

In Figure Bl the absolute deviation between the average 
formal velocity dispersion of an ensemble of 100 models and 
one with about 10000 models, which we take as a close to 
exact representation of the lens model, is shown in terms of 
percentage. The analysis is done for a subsample of doubles, 
quads and multiple object systems of SLAGS and GASTLES 
lenses. We finde that (TienB(-RionEi) (o"icna(-Rott)) of the smaller 
ensemble deviates less than ~ 1% (~ 2%) from the corre- 
sponding velocity dispersion for a larger ensemble. 



APPENDIX A: 

For a Hernquist profile the velocity dispersion (Ti, of a pro- 
jected distribution is 



aMi^^, [(2 + s^)X{s) - 3] i 2 (1 - s^)3 
X [~3s'^X{s){8s^ - 28/ + 35s^ - 20) 

^245^^ + 68.S* + 65s^ + 6] - 67rs I (Al) 



with s = R/a, whereas a denotes a scale length and 



© 2009 RAS, MNRAS 000, [THS] 



A lensing view on the fundamental plane 13 



Lens 




(J 


3bs 


0"lena 














Li 


Mie„./Ll 


Ref. 






[km 




[km 


[kpc] 


[kpc] 


[10^1 X Mq] 


[10" 


X 


Mq] 


[10" X Lq] 


[Mq/Lq] 




B0047-2808 


0.4850 


229 


± 15 


189.0 it 6.4 


5.324 


7.440 


2.084 ± 0.141 


3.059 


± 





401 


0.760 


2.74 


CD@® 


CFRS03.1077 


0.9380 


256 


± 19 


306.6 it 22.0 




15.905 
















CD®® 


Q0957 (2D) 


0.3600 


288 


± 9 


351.6tlll 


22.64 


25.529 


.^.^ + 2 477 
30.458if^^^ 


20.58 


-1- 


1 


29 


4.922 


6.19 


®®® 


PG1115+080 


0.3100 


281 


± 25 


190.8 ± 8.4 


2.072 


6.188 


0.826 ± 0.073 


1.792 


± 





319 


0.671 


1.23 


®(J)® 


HST14176 


0.8100 


230 


± 14 


245.6 ± 4.0 


5.190 


12.780 


3.430 ± 0.112 


3.008 


± 





366 


3.910 


0.88 


®®® 


HST15433 


0.4970 


108 


± 14 


156.2 ± 2.8 


- 


4.601 


- 




- 






- 


- 


®®® 


B1608+656 


0.6300 


247 


± 35 


242.6 ± 18.4 


4.291 


15.542 


2.766 ± 0.420 


2.868 


± 





813 


3.354 


0.82 


®®® 


MG2016+112 


1.0100 


304 


± 27 


330.4 ± 29.6 


1.707 


19.388 


2.041 ± 0.366 


1.729 


± 





307 


2.128 


0.96 


®®® 


Q2237+030 


0.0400 


215 


± 30 


146.4 ±2.8 


2.993 


0.743 


0.703 ± 0.027 


1.516 


± 





423 


0.542 


1.30 


®®® 


J0037-094 


0.1954 


265 


± 10 




6.804 


6.470 




5.234 


± 





395 


1.990 


1.99 


O 


J0737+321 


0.3223 


310 


± 15 


233.6 ± 6.4 


9.823 


5.333 


5.874 ± 0.322 


10.34 


± 


1 


00 


3.544 


1.66 


o 


J0912+002 


0.1642 


313 


± 12 


276.ot;i;= 


9.203 


5.271 


7 679+'' <'^'' 
' ■"'y_0.S21 


9.877 


± 





757 


2.289 


3.35 


o 


J0956+510 


0.2405 


299 


± 16 


266. St^;^ 


8.607 


6.201 


°' '-^"-0.414 


8.430 








2.572 


2.61 


o 




0.2150 


235 


i 10 


Tin O + 10.6 


7. 805 


5. 138 


. -,,+0.427 


4.722 


± 





402 


1.936 


2.34 


o 


J1330-014 


0.0808 


178 


± 9 


142.3t2«? 


1.244 


1.696 


0.276t«;°?^ 


0.432 


± 





044 


0.147 


1.88 


o 


J1636+470 


0.2282 


221 


± 15 


230.7t^2% 


5.256 


6.150 




2.812 


± 





382 


1.376 


2.23 


o 


J2300+002 


0.2285 


283 


± 18 


239.8t^"g'^ 


6.256 


5.278 




5.489 


± 





698 


1.522 


2.59 


o 


J2303+142 


0.1553 


260 


± 15 


242 7+l'l-9 


7.901 


5.303 


5 098+'' *""' 


5.851 


± 





675 


2.333 


2.19 


o 


AGO 1689 


0.1830 


1400 


± 300 


1 ISS 7+40.0 
1188.7_5g Q 


400.0 


237.6 


g^g2 + 424 


8589 


± 


2 


187 


33.48 


185 


e 


AGO 2667 


0.2330 


960 


+ 190 
-120 


762.0i^;^ 




98.01 


















Table Bl. 


Full set of 


gravitational lenses used 


for this 


analysis. 


The first 9 lenses are 


from CASTLES, the following 


9 from 



and the last two a re cl usters. The image po sitions and flux data have been taken from HST data (www.cfa.harvard.edu/glensdata/), 
iBolton et al.l | |2006| ) and lCovone et a h l |2006|). Th e sy mbols mark the references for t he data in the columns zi^, (Tabs and iJott and refer to 
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® llTonrvl 19981). ® llohvama et al.i.200a). (D iKoopmans et al.ll2003h. ® jlKoopmans fc Treul l2002'). ® (F oltz et al..i,1992 '). ® (.Keeton et al, 
Il998h . O Polton et al.l200d ). » l|Lokas et al.l200g) and ® l|C ovone et al.l2006r) . Note that (Tobs is the kinematic central velocity dispersion. 



which is in the case of SLACS lenses the line of sight stellar velocity dispersion measured by the 3" diameter SDSS spectroscopic fiber. 
The CTiona values are determined for the project e d dista nce Riena from th e oute rmost lensing image to the central lensing mass. The 
effective radii given in arcseconds in lRusin et al.l ll2003bl ) and lBolton et ahl ||2006|) have been transformed into kpc. All quantities in the 
table assume Ho = 72 km Mpc"'^, fljn = 0.3 and Qa = 0.7. 
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Table B2. List of previously found FP parameters and the results of this paper. 
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The following Figures |B2[ [B3] and [641 are for online publication only. 

Figure B2. CASTLES lenses used in this analysis. First column includes the Pixelens input files, the second column shows the formal 
velocity dispersion curves and in the third column the projected mass distributions (red dots mark the image positions, blue dots mark 
the source position) are presented. IMPORTANT NOTE: The y-axes of the velocity dispersion plots need to be multiplied by ^/^/ir ^ 0.8 
to yield the true aiens values. 



object B0047-2808 
symm pixrad 10 
redshifts 0.485 3.6 



1.270 
-0.630 

0.520 
-0.730 
g 13.7 



0.105 
-0.995 
-1.045 
0.705 




object CFRS03. 1077 
pixrad 10 

redshifts 0.938 2.941 

-2.085 0.383 
-1.979 -0.521 
-1.415 -1.213 
0.713 0.340 
g 13.7 




object Q0957-t-561 (2D) 
symm pixrad 10 
redshifts 0.356 1.41 
shear -30 
double 

1 . 408 5 . 034 
0.182 -1.018 423 
double 

2.860 3.470 
-1.540 -0.050 
g 13.7 




object PG1115-I-080 
pixrad 10 

redshifts 0.31 1.72 

0.381 1.344 
-0.947 -0.690 
-1.096 -0.232 
0.722 -0.617 
g 13.7 
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pixrad 10 
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g 13.7 
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pixrad 10 

redshifts 1.01 3.27 
double 

-1.735 1.778 
1.269 0.274 
g 13.7 



object (J2237-I-030 
pixrad 10 

redshifts 0.04 1.69 

0.598 0.758 
-0.075 -0.939 
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Figure B3. SLACS lenses used in this analysis. First column includes the Pixelens input files, the second column shows the formal 
velocity dispersion curves and in the third column the projected mass distributions (red dots mark the image positions, blue dots mark 
the source position) are presented. IMPORTANT NOTE: The y-axes of the velocity dispersion plots need to be multiplied by ■\/2/'k ^ 0.8 
to yield the true (Jicns values. 
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object J0737+321 
pixrad 10 
zlens 0.3223 
multi 3 0.5812 
-0.08 1.17 1 
0.52 -0.77 2 
-0.68 -0.57 2 
g 13.7 



object J0912 
pixrad 10 

redshifts 0.1642 0.3239 
double -1.275 1.025 0.375 
double -0.425 1.675 0.225 
double 0.975 1.575 -0.425 
double 1.475 1.225 -0.625 
g 13.7 




object J0956+510 
pixrad 10 
zlens 0.2405 
multi 3 0.47 
-1.67 0.170 1 
-0.220 -1.270 2 
0.730 0.620 2 
g 13.7 




object J1205+491 
pixrad 10 
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1.480 -0.320 
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object J2303+142 
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Figure B4. Lensing clusters used in this analysis. IMPORTANT NOTE: The y-axes of the velocity dispersion plots need to be multiplied 
by y/^/n ^ 0.8 to yield the true fJions values. 
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pixrad 18 
maprad 95 
minsteep 
shear 
g 13.7 
models 100 
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pixrad 10 
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models 100 



ACQ 1689 - continued - 
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